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摘 要 : 区 域 植被 动态 变化 研究 是 国家 级 自然 保护 区 生态 环境 效益 评估 的 重要 手段 之 一 。 基 于 谷 
歌 地 球 引 人 擎 (Google Earth Engine, GEE) Z F £ , 1E H Landsat 遥感 影像 提取 1995 一 2020 年 内 蒙古 大 
青山 国家 级 自然 保护 区 逐年 植被 生长 季 归 一 化 植被 指数 (Normalized difference vegetation index, 
NDVI) 数 据 集 ,采用 皮尔 还 相关 性 、 残 差分 析 和 Lindeman-Merenda-Gold(LMG ) 模 型 等 方法 对 该 保 
护 区 设立 前 后 的 NDVI 时 空 变 化 特征 及 其 影响 因子 进行 了 分 析 。 结 果 表 明 :(1) 1995—2008 年 NDVI 
呈 下 降 趋 势 的 区 域 面 积 点 比 为 69.04% , 而 在 2008 一 2020 年 NDVI 呈 上 升 趋势 的 区 域 占 比 高 达 
94.98% ,研究 区 成 为 国家 级 自然 保护 区 以 后 植被 明显 好 转 。(2) 1995 一 2008 年 负 向 人 类 活动 与 气 
候 暖 于 化 趋势 导致 实验 区 植被 NDVI 的 减少 ,其 中 气候 暖 干 化 趋势 是 主导 因素 (3) 降水 量 的 增加 
和 正 向 人 类 活动 驱动 了 2008 一 2020 年 实验 区 植被 NDVI 的 增加 。 其 中 ,以 生态 环境 保护 工程 为 代 
表 的 正 向 人 类 活动 是 该 时 间 段 植被 恢复 的 主要 原因 。 研 究 结果 可 为 大 青山 国家 级 自然 保护 区 环 
境 治 理 与 生态 工程 建设 提供 科学 依据 。 
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自然 保护 区 建设 是 我 国生 态 文 明 建设 的 重要 ”状况 的 重要 指标 ,能 够 定性 、 定 量 评价 植被 生长 和 
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组 成 部 分 ,在 维持 生物 多 样 性 Pi ZK A LE fie SE 
方面 起 着 关键 作用 "。 尤 其 是 近 些 年 来 ,我 国 在 国 
家 级 自然 保护 区 建设 和 管理 层面 上 取得 了 巨大 发 
展 ,截至 2018 年 ,已 建立 国家 级 自然 保护 区 474 处 ， 
对 区 域 目 然 资源 合理 利用 和 生态 环境 的 恢复 与 保 
护 方面 均 起 到 了 积极 作用 ”*。 由 此 ,自然 保护 区 生 
态 环 境 效 益 和 保护 成 效 评估 人 研究 也 成 为 诸多 学 者 
所 关注 的 热点 问题 之 一 一 。 

植被 动态 变化 研究 是 区 域 生态 环境 效益 评估 
的 有 效 手 段 。 归 一 化 植被 指数 (Normalized differ- 
ence vegetation index , NDVI1) 作 为 反映 地 表 植 被 生长 
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履 盖 状况 ,得 到 了 广泛 应 用 一 。 因 高 分 辨 率 遥 感 影 
像 使 用 门槛 较 高 ,已 有 研究 大 多 采用 中 、 低 分 辨 率 
影像 数据 开展 长 时 序 的 植被 NDVI 变 化 研究 ,但 
此 举 在 小 区 域 尺度 往往 无 法 满足 精度 上 的 需求 。 
Landsat 是 一 种 常见 的 中 分 状 率 遥感 影像 ,具有 时 间 
跨度 长 .多 波段 等 优势 ,在 免费 开放 的 影像 中 具有 
较 高 的 空间 分 状 率 ,在 植被 NDVI 获 取 方 面具 有 和 较 
高 的 实用 价值 。 然 而 , 因 该 影像 的 处 理 过 程 复杂 ， 
较 难 实现 长 时 间 序 列 的 植被 监测 工作 ”" 。 和 谷歌 地 
球 引 擎 (Google Earth Engine, GEE ) 云 平台 ,兼顾 了 
数据 收集 下 载 和 强大 的 空间 分 析 处 理 能 力 ,能 够 文 
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持 长 时 间 序 列 的 遥感 影像 获取 和 分 析 ,实现 了 数据 
获取 、 处 理 、 分 析 与 应 用 的 一 体 化 ,为 长 时 序 .中 高 
分 辩 率 植被 遥感 信息 提取 提供 了 新 的 技术 手段 "2 。 
阴山 山脉 是 中 国 北方 重要 的 自然 和 生态 界 山 ， 
其 中 大 青山 地 处 阴山 山脉 中 段 ,是 生态 多 样 性 最 丰 
富 的 区 域 之 一 。 该 保护 区 具有 涵养 水 源 、 水 土 保 
持 、 防 风 固沙 、 固 碳 释 氧 等 生态 功能 ,也 是 黄河 上 中 
游 重 要 的 水 源 补给 区 ,在 调节 我 国 北方 水 分 平衡 和 
水 资源 供给 中 起 着 重要 作用 ”。 然 而 ,气候 干旱 、 
土地 贫 交 ,加 之 人 为 过 度 开发 利用 导致 该 区 生态 环 
境 极 其 脆弱 …。 为 了 扭转 此 局 面 ,政府 部 门 于 2000 
年 在 内 蒙古 自治 区 成 立 了 自然 保护 区 ,在 2008 年 设 
立 了 国家 级 自然 保护 区 ,并 实施 了 一 系列 生态 环境 
保护 工程 。 但 目前 有 关 该 区 生态 环境 效益 评估 还 
未 得 到 有 效 开 展 。 因 此 ,本文 借助 CER 平台 ,采用 
Landsat 影像 对 比分 析 了 大 青山 国家 级 自然 保护 区 
建立 前 后 的 NDVI 时 空 演变 特征 ,并 定量 评估 了 和气 
候 因 子 和 人 类 活动 对 植被 变化 的 影响 。 研 究 成 果 
可 以 探 清 大 青山 国家 级 自然 保护 区 成 立 对 该 区 带 


(a) 保护 区 位 置 与 高 程 
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图 1 研究 区 概况 


Fig. 1 Overview of the study area 


2 数据 与 方法 


2.1 数据 来 源 

本 文 所 用 的 遥感 数据 是 Landsat-5 TM , Land- 
sat-7 ETM+ 和 Landsat-8 OLI 时 间 序 列 数据 。 选 择 
1995 一 2020 年 4 一 10 月 植被 生长 季 云 量 低 于 10% 的 
影像 (图 2) ,利用 GEE 平 台 对 数据 进行 格式 转换 、 数 
据 拼接 与 裁剪 ,并 计算 NDVI 值 '”"。 由 于 Landsat-8 
卫星 传 感 大 和 Landsat-5、7 卫星 传 感 需 之 间 有 明显 
差异 ,根据 Roy 等 提供 的 截 距 和 偏 移 调整 Land- 
sat-8 各 个 波段 的 反射 率 , 实 现 Landsat-s. 7 8 协同 ， 


乌 兰 察 布 市 人 (b) 保护 区 功能 区 分 布 À 


来 的 生态 效益 ,同时 可 为 保护 区 未 来 生态 环境 管理 
决策 提供 基础 信息 与 技术 文 返 。 


1 研究 区 概况 


Pye th Ka LL Bs Be H A RAP KMF E SRR 
川 平原 北部 的 呼和浩特 市 .包头 市 和 乌 兰 察 布 市 境 
P (109°47'~112°17'E, 40°34'~41°14'N) ,总 面积 约 
3888.70 km ,东西 长 度 为 217 km ,南北 平均 宽度 为 
18 km。 保 护 区 属 中 低 山 地 ,山区 海拔 为 1012~2319 m, 
北 坡 平缓 南 坡 陡峭 (图 la)。 该 区 属于 温带 半 干 旱 
大 陆 性 季风 气候 ,年 均 气 温 3~4 % ,年 降水 量 400~ 
500 mm。 保 护 区 北 坡 易 受 之 古 高 原 干 燥 气 流 的 影 
啊 ,气候 寒冷 .干燥 ; 南 坡 由 于 山地 的 阻挡 ,气候 温 
Be 湿润 。 该 保护 区 是 山地 森林 AA, .山地 草原 和 
重要 水 源 涵养 地 为 主 的 综合 性 山地 和 森林 类 保护 区 ” 。 
本 文采 用 的 保护 区 界线 以 国务 院 批 准 的 文件 为 准 .于 , 
保护 区 按照 不 同 功能 ,可 分 为 核心 区 缓冲 区 和 实 
验 区 (图 1b)55。 


并 以 最 大 值 合成 法 得 到 1995 一 2020 年 NDVI 年 值 数 
据 集 ,空间 分 辨 率 为 30 m"”。 该 数据 集 NDVI 最 大 
值 均 集中 于 7 月 ,最 大 值 在 0.89~1.00 之 间 ,年 际 变 
化 可 对 比 性 较 强 。 

本 文采 用 的 气象 数据 为 中 国 地 面 气温 降水 日 
值 0.5°*x0.5° 格 点 数据 集 (V2.0)。 选 取 人 研究 区 
1995 一 2020 年 4 一 10 月 生长 季 的 9 个 像 元 数据 , 利 
用 Python 对 数据 进行 反 距 离 权重 插值 ,并 通过 平 
均 求 和 获取 了 人 研究 区 气温 、 降 水 量 年 值 数 据 集 , 空 
间 分 辨 率 为 30 mo 

土地 利用 数据 源 自 武汉 大 学 杨 杰 `. 黄 昕 教授 基 
F GEE 利用 Landsat 数 据 提 取 的 1985—2020 年 中 
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图 2 Landsat 影 
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Fig. 2 Specific date of the Landsat images 


国土 地 覆盖 数据 集 ,空间 分 辨 率 为 30 m”。 本 文 提 
取 了 1995 一 2020 年 的 草地 和 林地 面积 ,用 于 表征 研 
究 区 生态 环境 工程 实施 成 效 。 

2.2 研究 方法 

2.2.1 Sen 民 斜 率 与 Mann-Kendall 非 参 数 检 验 Sen 
氏 和 斜率 是 非 参 估算 线性 趋势 的 算法 ” ,本 研究 用 于 
获取 1995 一 2008 年 和 2008 一 2020 年 2 个 时 间 段 植 
被 NDVI 的 变化 趋势 ,其 计算 公式 为 : 


i<j (1) 


B= Median J 
J 


式 中 :B 为 趋势 度 ; x A x, 分 别 为 某 一 像 元 在 ;年 和 7 
年 的 值 。 使 用 趋势 度 8 判 断 时 间 序 列 趋势 的 升降 ， 


当 B>0 时 ,表示 植被 NDVI 呈 增加 趋势 ; 当 B<0 时 , W 
本 文采 用 Mann-Kendall 方 法 在 0.05 置信 水 平 上 

对 NDVI 的 变化 趋势 进行 显著 性 检验 ”1。 

2.2.2 相关 性 分 析 本 文采 用 皮尔 逊 相关 系数 测定 

NDVI 与 气温 和 降水 量 之 间 的 相关 关系 汪汪 ,计算 公 

式 为 : 


———— (2) 
Seo $ iy) 


TU ry NAKA Rs n AEE I BOSE Ys xy HIR 
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分 析 的 2 个 变量 ;xy 分别 为 其 样本 值 ; x AUR 
子 的 平均 值 ; 7 为 NDVI 的 平均 值 。 本 文 同 样 在 0.05 
置信 水 平 上 对 相关 性 进行 显著 性 检验 汪 。 
2.2.3 RADA 本 文 假设 NDVI 仅 受 人 类 活动 和 
气候 条 件 影 响 ,并 利用 残 差分 析 法 评估 人 类 活动 对 
NDVI 的 影响 站。 计算 公式 为 : 
y=atbx,+cx, (3) 
o = NDVI — NDVI yay (4) 
式 中 :7 为 因 变 量 NDVI 值 ; x, 、x, 为 日 变量 ;4 为 第 
BX b c ATV ARM; o NIRE ; NDVI seul] NDVI gay 
分 别 为 实际 NDVI 值 和 预测 NDVI 值 。 残 差 趋势 为 
正 表 示人 类 活动 对 NDVI 起 促进 作用 ,为 负 表 示人 
类 活动 对 NDVI 起 抑制 作用 o 
2.2.4 Lindeman-Merenda-Gold 模型 相关 性 分 析 
和 残 差分 析 可 以 反映 气候 因子 与 人 类 活动 对 NDVI 
的 影响 ,但 无 法 量化 因子 对 NDVI 谈 化 的 贡献 率 ”。 
此 本 文通 过 多 元 线性 回归 模型 的 方差 分 解 量化 
各 变量 对 NDVI 变 化 的 相对 重要 性 。 并 选择 计算 最 


分 多 元 线性 回归 中 不 同 相 关 回 归 量 的 页 献 率 ,由 R 
语言 包 “relaimpo” 实 现 。 页 献 率 越 大 说 明 变 量 对 植 
被 变化 的 影响 越 大 ,反之 则 相反 ”。 选 取 的 气候 因 
子 包括 气温 和 降水 量 ; 人 类 活动 由 林地 和 草地 面积 
表征 , 即 人 工 造林 封山育林 .退耕 还 林 还 草 和 生态 
移民 等 生态 环境 工程 所 带 来 的 成 效 。 


3 结果 与 分 析 


3.1 NDVI 时 空 变化 特征 

3.1.1 年 际 变化 将 征 如 图 3 所 示 ,1995 一 2008 年 保 
护 区 植被 NDVI 斜 率 为 -0.002 a ', 显 著 性 检验 Z 值 
为 -0.77, 说 明 NDVI 整体 呈 非 显 车 下 降 趋 势 。 从 
NDVI 值 的 变化 看 出 ,2001 一 2003 年 NDVI 出 现 持续 
增加 , 且 增 速 较 快 ;2008 一 2020 年 研究 区 植被 NDVI 
呈 上 升 趋势 ,斜率 为 0.011 a ,2Z 值 为 2.20, 表 明 该 时 
间 段 内 NDVI 呈 显著 上 升 趋势 。 整 体 而 言 ,保护 区 
成 立 以 前 NDVI 呈 非 显 著 下 降 趋势 。2008 年 国家 级 
自然 保护 区 成 立 后 NDVI 呈 显著 上 升 趋势 ,说 明 国 


为 密集 的 Lindeman-Merenda-Gold(LMG) 模 型 2 | 


1995_2008 年 NDVI  ----- 1995 一 2008 年 Sen 式 斜率 
— 2008 一 2020 年 Sen 式 斜率 
0.62 
ee y= 0.002x+5.17 
0,52 7 
= Te a eee 
D 0.47 
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家 级 保护 区 的 成 立 促进 了 研究 区 植被 恢复 。 


一 e 一 2008 一 2020 年 NDVI 
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注 :NDVI 为 归 一 化 植被 指数 。 下 同 。 
图 3 NDVI 年 际 变化 特征 


Fig. 3 Characteristics of interannual variation of NDVI 


3.1.2 空间 变化 特征 1995 一 2008 年 保护 区 NDVI 
减少 区 域 共 占 总 面积 的 69.04% ,其 中 显著 减少 区 域 
占 7.78% , 主要 分 布 于 保护 区 东西 延伸 的 山脉 以 北 
地 区 (图 4a~b、 表 1)。NDVI 增 加 区 域 占 总 面积 的 
19.90% ,其 中 显著 增加 区 域 占 1.55% ,主要 分 布 于 保 


护 区 南部 沿线 地 区 。 说 明 在 国家 级 保护 区 成 立 前 
NDVI 变化 以 非 显 著 减 少 为 主导 ,呈现 “北部 减 小 、 
南部 增长 ”的 空间 变化 特征 。 分 区 方面 , 1995 一 
2008 年 核心 区 缓冲 区 和 实验 区 的 非 显 闭 减 少 区 域 
占 比 分 别 为 60.26% . 64.97% Fil 60.29% , 显著 减少 区 
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(a) 1995 一 2008 年 NDVI 变 化 趋势 
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图 例 NDVI 变 化 趋势 
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(c) 2008 一 2020 年 NDVI 变 化 趋势 
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(b) 1995 一 2008 年 NDVI 变 化 显著 性 检验 


ar 显著 性 
C P A DM 非 显 著 增加 sees 


稳定 不 变 


Bote BE 显著 减少 


(d) 2008 一 2020 年 NDVI 变 化 显著 性 检验 


图 例 显著 性 
Cae Moen Dee OO 稳定 不 变 
RE 


Bren BE ce 


图 4 NDVI 变 化 趋势 及 显著 性 检验 空间 分 布 
Fig.4 Spatial distributions of NDVI variation trend and significance test 


表 1 不 同时 期 NDVI 变 化 趋势 统计 


Tab.1 NDVI variation trend statistics in different time periods 


年 价 分 区 非 显 著 减 少 显著 减少 稳定 不 变 非 显 著 增加 显著 增加 
面积 /km Ek% ”面积 km 占 比 /% 面积 [km 占 比 /% 面积 [km h% 面积 km 占 比 /% 
1995—2008 核心 区 659.86 60.26 89.11 8.14 138.99 12.69 193.20 17.64 13.79 1.26 
缓冲 区 526.21 64.97 53.80 6.64 93.96 11.60 126.92 15.67 8.99 1.11 
实验 区 1196.10 60.29 159.49 8.04 197.22 9.94 393.49 19.83 37.57 1.89 
保护 区 2382.17 61.26 302.40 7.78 430.17 11.06 713.61 18.35 60.35 1.55 
2008—2020 ”核心 区 35.56 3.25 1.07 0.10 46.96 4.29 678.99 62.01 332.37 30.35 
缓冲 区 19.62 2.42 0.76 0.09 23.80 2.94 472.64 58.36 293.06 36.19 
实验 区 28.86 1.45 1.97 0.10 36.39 1.83 1046.07 52.73 870.58 43.88 
保护 区 84.04 2.16 3.80 0.10 107.15 2.76 2197.70 56.51 1496.01 38.47 
域 占 比 分 别 为 8.14% 、6.64% 和 8.04% ,说 明 3 个 分 区 ”于 其 他 功能 区 。 


植被 NDVI 的 变化 程度 较为 接近 。 

2008—2020 年 NDVI 减少 区 域 占 总 面积 的 
2.26% ,其 中 显著 减少 区 域 仅 占 0.10% ,分布 于 保护 
区 中 心地 带 。 而 NDVI 增 加 区 域 占 总 面积 的 
94.98% ,其 中 显著 增加 区 域 占 比 为 38.47% ,在 保护 
区 东 西部 地 区 较为 明显 (图 4c~d)。 说 明 国 家 级 自 
然 保 护 区 成 立 后 植被 NDVI 整 体 呈 现 出 向 好 的 特 
征 。 核 心 区 缓冲 区 和 实验 区 的 显著 增加 区 域 分 别 
占 其 总 面积 的 30.35% 、36.19% 和 43.88%。 说 明 在 


3.2 气候 因子 对 NDVI 的 影响 
3.2.1 和 气候 因子 的 年 际 变 化 将 征 1995 一 2008 年 保 
护 区 降水 量 的 变化 速率 为 -4.737 mm*a', 整体 呈 非 
显著 下 降 趋势 (图 5a); 气 温 则 以 0.053 Ca 的 速率 
呈 非 显著 上 升 趋势 (图 $b) ,说 明 该 时 期 保护 区 气候 
较为 干旱 ,不 利于 植被 生长 与 恢复 。 其 中 2001 一 
2003 年 降水 量 持 续 增 加 气温 持续 降低 、 气 候 相 对 
湿润 为 该 时 间 段 内 NDVI 的 增长 提供 了 有 利 条 件 。 
2008 一 2020 年 降水 量变 化 速率 为 3.059 mm-a', 


2008 年 目 然 保护 区 成 立 之 后 ,实验 区 的 植被 恢复 高 


整体 呈 非 显著 上 升 趋势 。 气 温 上 升 速率 由 0.053 °C .an 
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(a) 降水 量 
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图 5 气候 因子 年 际 变化 特征 


Fig. 5 Characteristics of interannual variation of climatic factors 


减缓 为 0.019 % .a ,说 明 气 候 暖 干 化 趋势 有 所 组 
解 。 尤 其 是 降水 量 的 增加 对 植被 生长 .恢复 起 到 促 
进 作 用 ,使 该 时 期 NDVI 增 加 。 

3.2.2 气候 因子 与 植被 NDVI 的 相关 关系 由 图 6 
和 表 2 可 知 ,2 个 时 间 段 内 保护 区 降水 量 与 NDVI 的 
相关 性 在 空间 上 有 所 不 同 。1995 一 2008 年 期 间 , 
NDVI 与 降水 量 的 相关 系数 以 正 值 为 主 ,其 中 非 显 
著 正 相关 区 域 占 79.06%。 显 著 正 相关 区 域 占 研究 
区 总 面积 的 18.54% ,主要 集中 于 保护 区 西部 和 东北 
部 ,在 中 部 地 区 也 有 零星 分 布 (图 6a~b) ,说 明 在 这 
些 地 区 降水 量 对 植被 的 影响 较 大 。 核 心 区 缓冲 区 
和 实验 区 的 显著 正 相 关 区 域 分 别 为 15.86%、16.90% 


FU 20.69% ,说明 降 水 量 对 实验 区 的 影响 大 于 核心 区 
和 缓冲 区 。 

保护 区 成 立 后 ,NDVI 与 降水 量 呈 非 显 著 正 相 
关 区 域 占 总 面积 的 72.08%。 显 著 正 相 关 区 域 占 比 
为 24.57% ,主要 集中 于 保护 区 西部 ,在 南部 沿线 和 
东北 部 地 区 也 有 体现 (图 6c~d) ,说 明 降水 量 对 ND- 
VI 的 影响 范围 在 扩大 。 并 且 , 这 些 地 区 植被 NDVI 
在 该 时 间 段 内 显著 增加 (图 6d) ,也 表明 了 降水 量 的 
增加 是 植被 恢复 的 主要 原因 之 一 。 核 心 区 .缓冲 区 
和 实验 区 的 显著 正 相 关 区 域 分 别 为 19.13% 19.58% 
和 29.61% ,说 明 降 水 量 对 实验 区 的 影响 仍 大 于 核心 
区 和 缓冲 区 。 
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#2 NDVI 与 降水 量 之 间 的 相关 性 统计 


Tab.2 Statistic for correlation between NDVI and precipitation 


年 从 分 区 显著 负 相 关 非 显 著 负 相关 显著 正 相关 非 显著 正 相 关 

面积 /km 占 比 /% 面积 /km 占 比 /% 面积 /km? 占 比 /% 面积 /km? 占 比 /% 

1995—2008 ”核心 区 0.14 0.01 28.07 2.56 173.66 15.86 893.08 81.56 
缓冲 区 0.37 0.05 18.62 2.30 136.84 16.90 654.05 80.76 

实验 区 0.37 0.02 45.86 2.31 410.48 20.69 1527.15 76.98 

保护 区 0.87 0.02 92.56 2.38 721.02 18.54 3074.25 79.06 
2008—2020 ”核心 区 0.13 0.01 45.09 4.12 209.50 19.13 840.23 76.74 
缓冲 区 0.16 0.02 30.37 3.75 158.57 19.58 620.78 76.65 

实验 区 0.34 0.02 54.32 2.74 587.35 29.61 1341.86 67.64 

保护 区 0.63 0.02 129.77 3.34 955.50 24.57 2802.80 72.08 

(a) 1995 一 2008 年 与 降水 量 相 关系 数 (b) 1995 一 2008 年 与 降水 量 相关 性 显著 性 检验 


N 


0 40km 
Ld 
图 例 相关 系数 图 例 显著 性 
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界线 -1.0-0.6-0.4-0.2 0.0 0.2 04 0.6 1.0 界线 E 非 显著 正 相关 非 显著 负 相关 


(c) 2008 一 2020 年 与 降水 量 相关 系数 (d) 2008 一 2020 年 与 降水 量 相关 性 显著 性 检验 


人 


40 km 

图 例 相关 系数 图 例 显著 性 
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界线 -1.0-0.6-0.4-0.2 0.0 0.2 0.4 0.6 1.0 界线 FU 非 显著 下 相关 PERS HA 


图 6 NDVI 与 降水 量 之 间 的 相关 系数 及 显著 性 检验 空间 分 布 


Fig. 6 Spatial distributions of correlation coefficient and significance test between NDVI and precipitation 


如 图 7 和 表 3 所 示 ,1995 一 2008 年 保护 区 气温 ”的 影响 最 小 。 
与 NDVI 以 负 相 关 为 主 ,其 中 非 显 著 负 相关 区 域 占 3.3 人 类 活动 对 植被 NDVI 的 影响 
研究 区 总 面积 的 56.17%。 显 著 负 相关 区 域 占 比 为 根据 国家 相关 规定 ,保护 区 内 的 核心 区 和 缓冲 
41.75%, 在 研究 区 东北 、 中 部 和 西部 地 区 均 有 体现 。 区 禁止 人 类 活动 。 因 此 ,本 文 以 实验 区 为 例 , 评 估 
(图 7a~b) ,说 明 气 温 的 升 高 导致 这 些 区 域 植被 NDVI ”人 类 活动 对 NDVI 的 影响 。 图 8 和 表 4 所 示 ,1995 一 
的 减少 。 分 区 方面 ,核心 区 缓冲 区 和 实验 区 的 显 。” 2008 年 残 差 趋势 为 抑制 NDVI 的 区 域 占 实验 区 总 面 
著 负 相关 区 域 占 比 分 别 为 41.01% 、41.47% 和 积 的 33.84% , 主要 分 布 于 北部 沿线 地 区 , 均 为 NDVI 
42.27%, 说 明 气 温 对 3 个 功能 区 的 影响 较为 一 致 。 下 降 或 显著 下 降 地 区 ,说明 不 合理 的 人 类 活动 导致 

2008 一 2020 年 气温 与 NDVI 以 非 显 著 正 相关 为 。 实验 区 北部 NDVI 的 下 降 ; 残 差 趋势 为 促进 NDVI 的 
主 , 占 保护 区 的 84.46%, 显 著 正 相 关 区 域 仅 占 1.58% ”区 域 占 实验 区 总 面积 的 51.57% , 主要 分 布 在 中 部 与 
(图 7c~d) ,说 明 该 时 间 段 内 气温 对 NDVI 无 显著 影 。” 南部 (图 8a) ,与 NDVI 呈 增加 趋势 的 地 区 分 布 较为 
响 。 核 心 区 、 绥 冲 区 和 实验 区 显著 正 相 关 区 域 分 别 。 一 致 (图 4a) ,说 明正 向 的 人 类 活动 促进 了 该 区 植被 
FJ 2.40% 2.15% 、0.89% ,表明 气温 对 实验 区 NDVI ”NDVI 的 增加 。 
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(a) 1995 2008 年 与 气温 相关 系数 (b) 1995 一 2008 年 与 气温 相关 性 显著 性 检验 


图 例 相关 系数 图 例 显著 性 
| 保护 区 界线 DE | [保护 区 界线 EE 显著 正 相 关 ME 显著 负 相 关 
1.0 —0.6 -0.4 -0.2 0.0 02 0.4 0.6 1.0 P SEEK [| 非 显 著 负 相关 
(c) 2008 一 2020 年 与 气温 相关 系数 (d) 2008 一 2020 年 与 气温 相关 性 显著 性 检验 


关系 图 例 显著 性 
[保护 区 界线 M [保护 区 界线 BE rrn BE 显著 负 相 关 
1.0 —0.6 —0.4 —0.2 0.0 0.2 04 06 1.0 E 非 显 著 正 相关 [EBS MX 


图 7 NDVI 与 气温 之 间 的 相关 系数 及 显著 性 检验 空间 分 布 


Fig.7 Spatial distributions of correlation coefficient and significance test between NDVI and temperature 


表 3 NDVI 与 气温 之 间 的 相关 性 统计 


Tab.3 Statistic for correlation between NDVI and temperature 


Ep 分 区 显著 负 相 关 非 显 著 负 相关 显著 正 相关 非 显著 正 相关 
面积 /km 占 比 /% 面积 /km 占 比 /% 面积 /km 占 比 /% 面积 /km 占 比 /% 
1995—2008 ”核心 区 449.04 41.01 623.33 56.93 0.12 0.01 22.46 2.05 
缓冲 区 335.83 41.47 459.06 56.68 0.24 0.03 14.76 1.82 
实验 区 838.55 42.27 1101.95 55.55 0.22 0.01 43.15 2.17 
保护 区 1623.42 41.75 2184.33 56.17 0.58 0.01 80.37 2.07 
2008—2020 ”核心 区 0.08 0.01 134.77 12.31 26.26 2.40 933.84 85.29 
缓冲 区 0.11 0.01 98.12 12.12 17.38 2.15 694.27 85.72 
实验 区 0.38 0.02 309.42 15.60 17.67 0.89 1656.41 83.49 
保护 区 0.57 0.01 542.34 13.95 61.30 1.58 3284.50 84.46 
(a) 1995 一 2008 年 (b) 2008 一 2020 年 
À À 
0 40km 0 40km 
a} Lcd 


图 例 | | 实验 区 。” 残 差 变 化 趋势 ae ce yy) et 


图 8 实验 区 NDVI 残 差 趋势 空间 分 布 


Fig. 8 Spatial distributions of NDVI residual trends in the experimental area 


2008 一 2020 年 实验 区 残 差 趋势 为 抑制 NDVI 的  NDVI 呈 增加 趋势 的 区 域 基本 重 麦 (图 8b)。 说 明 
区 域 占 比 仅 为 2.79%; 而 残 差 趋势 为 促进 NDVI 的 区 2008 年 国家 级 自然 保护 区 成 立 以 后 人 类 活动 对 ND- 
域 明显 增加 ,和 履 盖 实验 区 94.62% 的 面积 ,与 图 4b 中 VI 的 影响 以 促进 作用 为 主 ,并 驱动 了 对 实验 区 植被 
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表 4 实验 区 NDVI 残 差 趋势 统计 
Tab. 4 Statistics of NDVI residual trends in the 


experimental area 


Jan 1995—2008 4 a 年 
面积 /km 占 比 /% 面积 /km? 占 比 /% 
抑制 671.34 33.84 55.35 2.79 
无 影响 289.46 14.59 51.40 2.59 
促进 1023.07 51.57 1877.12 94.62 
的 恢复 。 


3.4 气候 因子 和 人 类 活动 对 NDVI 变 化 的 贡献 率 
为 了 进一步 量化 气候 因子 和 人 类 活动 对 NDVI 
变化 的 影响 ,同时 也 为 了 找 出 保护 区 成 立 前 后 植被 


(a) 1995 一 2008 年 
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NDVI 变 化 的 主导 因素 ,本 文 利 用 LMG 模 型 计算 了 
气温 、 降 水 量 .林地 面积 和 草地 面积 对 NDVI 的 贡献 
率 。 结 果 显 示 ,1995 一 2008 年 降水 量 对 实验 区 NDVI 
变化 的 贡献 率 为 32.00%, 气 温 的 贡献 率 为 58.42%， 
林地 和 草地 面积 的 贡献 率 不 足 10%( 图 9a) ,说 明 气 
候 暖 干 化 趋势 是 该 时 间 段 实验 区 NDVI 减 少 的 主要 
影响 因素 ;2008 一 2020 年 林地 面积 对 实验 区 NDVI 
变化 的 贡献 率 为 34.28% ,草地 面积 的 贡献 率 为 
43.93% ,降水 量 和 气温 对 NDVI 的 贡献 率 分 别 为 
19.02% 和 2.78%( 图 9b)。 说明 大 青山 国家 级 自然 保 
护 区 成 立 后 , 正 向 的 人 类 活动 为 实验 区 NDVI 变 化 
的 主要 影响 因素 ,也 说 明 实验 区 内 实施 的 生态 环境 


(b) 2008 一 2020 年 
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图 9 气候 因子 和 人 类 活动 对 实验 区 NDVI 变 化 的 贡献 率 


Fig.9 Contribution of climate factors and human activities to NDVI change in the experimental area 


保护 工程 取得 了 显著 效益 。 
4 讨论 


本 文 借助 GCEE 云 平台 ,利用 1995 一 2020 年 Landsat 
遥感 影像 ,提取 了 研究 区 30 m 分 状 率 植被 NDVI 数 
据 集 ,对 大 青山 国家 级 自然 保护 区 长 时 序 植被 动态 
变化 进行 了 监测 。 结 果 显 示 在 国家 级 自然 保护 区 
成 立 后 ,NDVI 呈 增加 趋势 的 地 区 占 人 研究 区 总 面积 
的 94.98%, 植 被 生态 环境 大 幅 好 转 , 这 与 胡 尔 查 等 
和 王 佳 新 等 "的 研究 结果 较为 一 致 。 因 选用 了 
Landsat 遥感 影像 ,上 是 借助 了 GEE 集 成 运算 环境 ,本 
文 获取 了 时 间 跨 度 长 .空间 分 辨 率 较 高 的 植被 监测 
结果 ,更 为 精确 地 刻画 了 植被 NDVI 的 时 空 演变 特 
征 ,丰富 了 长 时 序 、 小 区 域 尺度 植被 监测 技术 手段 。 

本 文 研究 结果 显示 气候 的 暖 干 化 是 国家 级 自 
然 保 护 区 成 立 前 NDVI 减 少 的 主要 原因 ,这 与 陈 淑 
君 等 ”的 研究 结果 较为 一 致 。 气 温 升 高 、 降 水 减少 
容易 导致 气象 干旱 的 发 生 , 从 而 使 得 大 气 .土壤 水 


分 亏 缺 ,加 剧 植被 的 干旱 胁迫 ,抑制 植被 生长 ?。 
本 文 还 发 现 大 青山 国家 级 自然 保护 区 成 立 后 , 正 问 
的 人 类 活动 是 实验 区 NDVI 变化 的 主要 影响 因素 。 
据 不 完全 统计 ,保护 区 已 累计 完成 人 工 造林 6.45 x 
10° hm ,完成 封 育 面积 1.32x105 he, 这些 工程 措施 
使 得 大 青山 国家 级 自然 保护 区 和 森林 履 盖 率 大 面积 
提高 后 。 另 外 ,以 往 的 研究 大 多 采用 土地 利用 类 型 
转移 面积 ,如 耕地 转化 为 林地 和 草地 面积 ” ”表征 
生态 工程 实施 成 效 。 而 本 文 直接 采用 草地 和 林地 
面积 予以 表征 ,原因 在 于 研究 区 耕地 面积 较 少 ” ， 
耕地 转 为 林 草 地 面积 仅 能 表征 小 部 分 生态 环境 工 
程 效 益 。 

本 研究 聚焦 于 国家 级 自然 保护 区 设立 对 生态 
环境 产生 的 正 向 作用 ,研究 结果 也 证 明了 生态 工程 
实施 后 植被 问好 发 展 的 事实 。 但 是 由 于 缺乏 放牧 
强度 砍伐 树木 等 表征 负 向 人 类 活动 的 时 空 数据 ， 
本 人 研究 未 能 将 这 些 指 标 应 用 到 LMG 模型 的 运算 。 
也 因 和 缺乏 保护 区 生态 环境 工程 实施 的 准确 信息 ,未 
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紧密 结合 生态 工程 实施 的 具体 时 间 、 方 位 等 开展 
深层 次 的 人 类 活动 对 植被 NDVI 变 化 的 驱动 机 制 研 
究 。 此 外 ,本 文 只 讨论 了 气温 和 降水 量 两 种 气候 因 
子 对 NDVI 的 有 影响。 内 壹 古 大 青山 国家 级 自然 保护 
区 作为 半 干 旱地 区 ,植被 对 气候 因子 的 响应 关系 是 
一 个 复杂 的 过 程 ,植被 变化 还 会 受到 潜在 蒸 散 等 其 
他 气候 因子 的 影响 中。 因此 ,在 未 来 将 考虑 多 种 气 
候 因 子 对 植被 NDVI 的 综合 影响 。 


II》 


5 结论 


(1) 1995—2008 年 保护 区 成 立 之 前 ,植被 NDVI 
以 -0.002 ai 的 速率 减少 ,减少 区 域 共 占 研究 区 总 面 
积 的 69.04%。2008 一 2020 年 保护 区 成 立 以 后 植被 
NDVIVA 0.011 a 的 速率 增加 ,增加 区 域 占 总 面积 的 
94.98% ,说 明 保护 区 的 建立 对 植被 的 好 转 起 到 明显 
促进 作用 。 

(2) 1995 一 2008 年 研究 区 气候 呈 暧 干 化 趋势 。 
NDVI 与 气温 呈 显 著 负 相关 区 域 面积 占 比 为 
41.75% , 与 降水 量 呈 显著 正 相 关 区 域 面积 占 比 为 
18.54%。2008 一 2020 年 NDVI 与 降水 量 、 气 温 均 以 
正 相 关 为 主 ,显著 正 相 关 区 域 分 别 为 24.57% 和 
1.58% ,说 明 在 该 时 间 段 内 降水 量 对 植被 NDVI 的 影 
响 有 所 增加 ,而 气温 的 影响 则 显著 减弱 。 

(3) 1995 一 2008 年 降水 量 和 气温 对 实验 区 ND- 
VI 变化 的 贡献 率 分 别 为 32.00% 和 58.42% ,说 明 气 
修 因 子 是 实验 区 NDVI 变化 的 主要 影响 因素 。 
2008 一 2020 年 实验 区 林地 面积 和 草地 面积 对 NDVI 
变化 的 贡献 率 分 别 为 34.28% 和 43.93%, 说 明正 向 的 
人 类 的 活动 是 植被 恢复 的 主导 因素 。 
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Abstract: In recent years, China has made great progress in the construction and management of national nature 
reserves such as the Daqingshan National Nature Reserve in Inner Mongolia. However, the associated ecological 
benefits have not been effectively assessed. The purpose of this study is to assess whether the establishment of 
the Daqingshan National Nature Reserve has contributed to the ecological recovery and improvement of the area. 
Examining changes in vegetation dynamics can be an effective tool for regional ecological engineering assess- 
ment, and this study analyzes such changes using the normalized difference vegetation index (NDVI). Based on 
the Google Earth Engine (GEE) cloud platform, the study uses Landsat remote sensing imagery to extract NDVI 
data for the nature reserve for the years 1995 to 2020. Spatial and temporal variation in NDVI and variation driv- 
ers before and after the establishment of the reserve were analyzed using Pearson correlation, residual analysis 
and the Lindeman-Merenda-Gold (LMG) model. A decreasing NDVI trend was identified during the 1995—2008 
period in 69.04% of the studied area, while an increasing NDVI trend was identified during the 2008—2020 peri- 
od in 94.98% of the studied area. These results indicate that the quality of vegetation in the study area has im- 
proved significantly since it became a national nature reserve. Negative impacts from human activities and climat- 
ic warming during the 1995—2008 period led to decreases in NDVI in the studied area, with climate warming be- 
ing the dominant factor. An increase in precipitation and positive impacts from human activities drove increases 
in NDVI in the studied area during the 2008—2020 period, with positive human impacts arising from ecological- 
environmental protection engineering implementation being the main reasons for vegetation recovery during this 
period. The selection of Landsat remote sensing images and the use of the GEE integrated computing environ- 
ment enabled the study to obtain vegetation monitoring data over a long time span and at a high spatial resolu- 
tion. Spatiotemporal variation in vegetation was also more accurately portrayed using the NDVI measure, thereby 
enriching the technical means for long time-sequence and small region-scale vegetation monitoring. The results 
of this study provide clear evidence of the ecological benefits brought by the establishment of the Daqingshan Na- 
tional Nature Reserve. At the same time, the study provides basic information and technical support for future 
ecological-environmental management decisions concerning the protected area. 
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